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ABSTRACT 



A laboratory study has been made of the characteristics 
of an interf erometric type optical fiber pressure gradient 
hydrophone. The optical system is configured as an 
all-fiber Mach-Zehnder interf erometer excited by a single 
mode gallium arsenide diode laser. A pair of identical 
fiber coils, one in each arm of the interferometer , forms 
the sensing portion of the gradient hydrophone. Each coil 
consists of 10 meters of polyethylene jacketed, single mode 
fiber, wound in a doughnut shaped element of mean diameter 4 
cm and thickness 3 mm. An acoustic calibrator, similar in 
design to the USRD type 640 shipboard calibrator used by the 
U.S. Navy, was constructed for use in obtaining single and 
dual (dipole) coil sensitivities. The acoustic sensitivity 
of each coil was determined separately and then the pair was 
aligned coaxially and separated by 10 cm to operate as a 
pressure gradient device. Details of the construction of 
the system, calibration procedure, lock-in demodulator, and 
hydrophone sensitivity data are presented. The system has 
demonstrated a 14 dB sensitivity advantage over the USRD 
type 663 standard (piezoelectric) pressure gradient 
hydrophone. 
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I 



INTRODUCTION 



A. BACKGROUND 

Th* usa o-F optical fibarB for communication and sanaing 
is a rapidly axpanding araa of racant tachnology. Optical 
fibar tachnology is a little ovar a dacada old, starting 
with tha first damonstrations of low-loss (lass than 1000 
dB/kffi) optical fibar wavaguidas. Tachnological advancas in 
matarial procassing, fabrication of optical fibars, 
intarconnaction davicas, cablas, sourcas, and datactors now 
allow repaatarless light transmission ovar graat distancas 
<tans of kilomatars) with an axtramaly larga bandwidth CRaf. 

n. 

Capitalizing on tha availability of optical fibar 
components, there has been significant interest and progress 
in the past half dacada in exploitation of optical fibers as 
transduction al aments in a new class of sensor s. These 
sensors are capable of detecting acoustic fields, linear and 
rotational acceleration, temperature, strain, electric and 
magnetic fields, and other physical parameters. This 
versatility can also lead to difficulties and much effort 
has baen expended to enhance sensitivity to the parameter of 
interest while desensitizing the fibar to all other 
fluctuations CRef. 21. Fiber sensors offer a number of 
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geometric 



potential advantages: increased sensitivity, 

versatility as -fiber sensors can be con-figured in arbitrary 
shapes, lightweight components, and low cost. Progress in 
demonstrating these advantages has been substantial in the 
past -few years with over 60 di-fferent types o-f sensors using 
optical -fiber as the transduction element being developed. 
These sensors range -from simple on /o-f -f fluid level 
indicators, to the more sophisticated interf erometrlc 

configurations. This large number of individual devices is 
usually categorized into amplitude or phase sensors. In 
effect, the sensor modulates some parameter of the lightwave 
in an optical fiber, such as the intensity or the phase. In 
the former case the physical perturbation interacts with the 
fiber to directly modulate the intensity of the light in the 
fiber. In the latter case the physical perturbation 
modulates the optical phase of coherent light propagating 
through the fiber. Optical interf erometry must be utilized 
to convert the optical phase modulation to optical intensity 
modulation in order to permit photodetection. This 

technique allows for extremely high sensitivity. Phase 

modulation, conversion of phase modulation to intensity by 
interf erometry , and subsequent detection by photodiodes will 
be discussed in more detail in Chapter II. 
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B. PURPOSE OF THIS STUDY 



Fib»r optic acoustic ssnsors ware -First dsmonstrated to 
bs -feasibla in 1977 CRa-F. 3 and 41. Sinca than a signi-Ficant 
dagraa o-F programs has baan achiavad, both in tha 
undarstanding o-F tha ralavant acousto-optic transduction 
machanisms as wall as in componant davalopmant and sansor 
packaging CRa-F. 9 and 61. Most of tha fibar optic hydrophona 
work par for mad to data has baan with monopolar 
omni dir actional davicas. Figura 1.1 illustratas tha basic 
gaomatry of a typical fibar optic hydrophona which utilizas 
acoustically inducad phasa modulation for transduction. A 
lasar baam is split into two arms of an intarfaromatar whara 
a sansing coil in ona arm is axposad to tha acoustic 
prassura fiald whila a rafaranca coil, forming tha sacond 
arm, is i sol at ad from tha fiald. Tha two baams ara 
racombinad and allowad to intarfara on tha surfaca of a 
photodatector . 

Sinca thair invantion in tha lata 1970‘’s two unique 
char actari sties of intarf aromatric typa fibar optic 
hydrophonas, i.a. , thair basic two-arm dasign and thair 
provan high sansitivity, hava suggastad thair usa as 
gradiant-typa davicas. A fibar optic gradiant hydrophona 
has a sensing coil in each arm of tha intarf aromatar and 
both are exposed to the pressure fiald. Figura 1.2 
illustrates the basic geometry of an acoustic fiber optic 
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Figure 1.1 Fiber Optic Interferometric 
Hydrophone 
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SENSOR 



• BOTH SENSOR AND REFERENCE COIL 
INSONIFIED 

• DETECTS SPATIAL VARIATION IN 
EXCITATION FIELD 



Figure 1.2 Fiber Optic Interferometric 
Gradient Hydrophone 
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gradient hydrophone. It was the aim of this project to 
construct and partially evaluate in a laboratory environment 
the character! sties of a fiber optic gradient hydrophone. 

C. FORMAT OF THE REPORT 

The theoretical bases of conventional gradient 
hydrophones, the calibration of gradient hydrophones, and 
the behavior of interf erometrlc sensors used in preparation 
of this report is presented in Chapter II. Details of the 
construction of the optical fiber sensor are presented in 
Chapter III. Included in Chapter III is a description of the 
characteristics of the individual parts of the 
interf erometer . The parts are discussed in the order that 
they were built into the system. Specific instrumentation 
used to obtain data is also discussed in Chapter III. The 
experimental procedures used to establish the 
characteristics of the system and conduct sensitivity 
measurements of the individual sensor colls and the gradient 
(dual coil) hydrophone are discussed in Chapter IV. Analysis 
of the data and evaluation of the results are also presented 
in Chapter IV. Chapter V contains concluding remarks and 
recommendations for further work. Appendix A lists a BASIC 
computer language program used in data acquisition and 
Appendix B presents an example of the data obtained using 
this program. 
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I I . THEORY 



A. CONVENTIONAL GRADIENT HYDROPHONE 

In many acoustic applications it is necessary to 
determine the direction o-f signal arrival as well as to 
detect acoustic pressure levels. Directional -features are 
usually achieved with the use o-f multiple hydrophones 
spatially distributed in some speci-fied orientation, ■•g«, a 
line array. The simplest o-f these designs is the dipole 
pressure gradient hydrophone. 

Some pressure gradient hydrophones are designed so that 
the sensing part o-f the hydrophone oscillates in the medium 
under the in-fluence o-f the sound pressure gradient. That 
is, a hydrophone is constructed -for which the driving -force 
and the hydrophone's output is proportional to the 
difference between the pressures at two di-f-ferent locations 
in an acoustic sound -field. Ther e-fore, the electrical 
output is proportional to the pressure gradient in the sound 
-field. The main characteristic and advantage o-f a pressure 
gradient hydrophone is the directivity pattern obtained -from 
such a small .size hydrophone. As will be shown below, 
pressure gradient hydrophones have dipole, or figure-eight, 
directivity patterns, as indicated in Figure 2.1; hence they 
are bidirectional. Assuming the size of the hydrophone is 
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0» 




Figure 2.1 Directivity Pattern of Pressure 
Gradient Hydrophone 
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9 



the dipole 



small compared to the acoustic wavelength \ 
response when oriented at any angle Q relative to an 
incoming plane pressure wave is proportional to cos ^ . 

Both for simplicity and since they resemble the fiber 
optic system under study, pressure gradient hydrophones 
formed by combining two small, closely spaced, individual, 
pressure hydrophones will be discussed. The performance of 
this hydrophone pair (or doublet) is easily calculated, and 
the behavior thus derived applies also to pressure gradient 
hydrophones of more complicated construction as long as the 
assumptions of compactness made in the derivation are not 
violated. 

As an introduction to the theory of this type of 
hydrophone, consider two small pressure hydrophones placed 
in a standing acoustic wave field, P(x,t), separated by a 
distance d << \ as shown in Figure 2.2. The hydrophones are 
placed so that a null in the standing wave pattern occurs 
halfway between the two hydrophones. The electrical outputs 
of the two are connected 180 degrees out-of— phase with each 
other to obtain a dipole. The dimensions of the two 
hydrophones are typically much less than a wavelength, 
therefore the acoustic field is assumed to be negligibly 
disturbed. 
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M|q. 




0 ^ angle of incidence 

d — spacing between hydrophones 

M — Individual hydrophone free 
field voltage sensitivity 



Figure 2.2 Geometry Used in Deriving Sensing 
Characteristics of Acoustic Dipole 
(Pressure Gradient) 
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Consider the plane sinusoidal wave shown in Figure 2.3. 
The instantaneous acoustic pressure, P(x,t) is given by: 

jCJt 

P(x,t) « P BinCkxle 
o 

( 2 . 1 ) 

where P is the peak acoustic pressure, k is the propagation 
o . 

wave number k = 2 7T/Af ^nd x is the distance o-f one o-f the 
hydrophones -from absolute minimum pressure (i.e. the 
standing wave pressure node). Assuming the argument of the 
sine function is small, the first order (linear) 
approximation (i.e., sin y * y) of the pressure amplitude at 
X is given by: 



jOJt 

P(x,t) - P kx e 
o 

( 2 . 2 ) 

If one sensor is at x « -d/2 and the other at x » as 

shown in Figure 2.3, then the amplitude of the pressure 
difference to be measured may be expressed as: 



Ap = P - P 

+d/2 -d/2 



(2.3) 
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Figure 2.3 




Pressure Distribution as Function of 
Distance from Null Pressure Point 



22 



This may be written as: 



Ap 



P kCd/2 - (-d/2) : 



o 



(2.4) 



Or: 



Ap 



P kd 



o 



(2.5) 



For the more general one dimensional case o-f any point 
along a standing wave, consider a standing wave in a tube o-f 
water as shown in Figure 2.4. The instantaneous acoustic 
pressure P(x,t) is given by: 



where x has been chosen to be zero at the air-water 
boundary. 

If we expand equation (2.6) into a Taylor's series about 
X = h, where h is the distance from the air— water boundary, 
the following equation results: 



JCJt 



P(x,t) = P sin(kx)e 



o 



( 2 . 6 ) 




2 



h) + 



(2.7) 
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h=0 



X 



± 



Figure 2.4 Pressure and Velocity Distribution in a 
Standing Wave Tube. 
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where the -first term and all their partial derivatives are 
evaluated at x sh. Timet is held constant. 

The -first partial derivative (pressure gradient) is: 



dp/d 



X 



jCJt 

P kcos(kx)e 
o 



( 2 . 8 ) 



The second partial derivative is: 

^2^2 2 jGJt 

U P/C/x » -P k sin(kx)e 
o 

( 2 . 9 ) 

The pressure di-f-ference between two points axially 
displaced -from one another by distance ^x where ^x = x - h 
is: 

Ap = P(x,t) - P(h,t) 



( 2 . 10 ) 



or 



Ap = <3 p/3 X) Ax + 



( d p/d X ) (Ax) 



2 

/2 + . . . 



( 2 . 11 ) 

I-f Ax is small compared to the wavelength o-f the standing 
wave all terms above second order may be neglected and 
there-fore the series may be truncated a-fter the second order 
term. 
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I-f h and ara salactsd so that a prassura noda is 
locatad at Ax/2, than tha diffaranca in prassura Ap 
raducas to aquation (2.S> with Ax > d, in tha linaar 
approximation. Tha first corraction, 6 , to tha linaar 
approximation, in tha vicinity of a prassura noda, bacomasi 

6 " -C <kAx)/23tan(kh) 



( 2 . 12 ) 

A maximum in tha prassura diffaranca will occur whan 
aach sansor is aqually distant from tha prassura noda and 
convarsaly a minimum in tha prassura diffaranca will occur 
whan aach sansor is aqually distant from a prassura maximum 
(antinoda) . 

A prassura grad i ant hydrophona is somatimas cal lad a 
velocity hydrophona, the implication being that it responds 
to the particle velocity in tha acoustic medium. This is 
somewhat misleading since the velocity of tha hydrophone 
surface is normally much smaller than the particle velocity 
of the incident acoustic wave. But even though tha 
hydrophone is not a direct detector of particle velocity, 
its response, being proportional to prassura gradient, is 
also proportional to particle velocity of tha incident 
acoustic wave as shown in the one dimensional linearized 
Euler aquation: 

-0P/0X - jCjpu (2.13) 
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where CJ is the angular frequency of the acoustic wave, p 
the mass density of the medium, and u is the Eulerian 
particle velocity of the free-field acoustic signal. 

B. CALIBRATION OF ORADIENT HYDROPHONES 

Pressure gradient hydrophones may be, and usually are, 
calibrated in terms of pressure rather than velocity 
(pressure gradient is never used as a reference). Since the 
relationship between pressure and velocity depends on the 
wavs and boundary conditions, these conditions must be 
known, be standard, or fit the definition of sensitivity. 
Therefore the sensitivity of pressure gradient hydrophones 
usually is given in terms of volts/micropascals (V/pPa) 
specified at a particular frequency. 

Plane progressive waves are specified in the definition 
of free-field voltage sensitivity, and the relationship 
between the pressure p and velocity u in such waves is given 
by the wave or characteristic impedance p/u = pc. Thus, the 
sensitivity of a pressure gradient hydrophone expressed in 
terms of pressure differs from the sensitivity expressed in 
terms of velocity by the constant factor pc, where p is the 
density of the medium and c is the speed of sound. 

The standard procedure for pressure gradient hydrophone 
calibration is by direct comparison with a standard pressure 
hydrophone only if essentially plane waves impinge on both 
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hydrophones. Since waves in a free field never are 
perfectly plane or if a spherically diverging wave from a 
point source is used, a correction factor. 



.2 2 

lOlogCl-t- A /(27rr) 1, where r is the distance from the 
source, must be subtracted from the measured sensitivity 
CRef. 71. This correction factor decreases as r increases. 

Because of the difficulties in obtaining free-field 
conditions particularly at low frequenclas, a standing-wave 
tube like that shown in Figure 2.4 was employed for the 
laboratory investigations described in this report. Since a 
water-air boundary is a near zero impedance boundary, 
essentially complete reflection will take place at the top 
of the tube. A standing-wave system then can be established 
and the following relations hold in the ideal case (i.e. 



SWR- OO ) I 



p ■ p sin kh, 
o 



(2. 14) 



u = (p /pc) com kh, 
o ^ 

(2. IS) 



p/u ^ pc tan kh. 



(2. 16) 
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If the pr«BSur« is measured at a point a distance h from 
the water— air surface, equations (2. 14) and (2. 15) can be 
used to compute the pressure and velocity at any other 
point. As stated earlier the hydrophones are assumed to 
have a negligible effect on the standing wave pattern. If 
the tube is operated below cutoff for non-ax ial modes, then 
the pressure variations normal to the tube axis can be 
neglected. A detailed discussion of the calibrator used is 
given in Chapter III section D7. 

A more detailed explanation of sensitivity and 
calibration procedures can be found in Bobber CRef 71. In 
this experiment sensitivities of single coil fiber optic 
acoustic hydrophones are defined in terms of ^ radians/^U Pa. 
Sensitivities of gradient (dual coil) hydrophones are 
defined in terms of optical phase shifts induced by acoustic 
pressures and are expressed in /i radians/^ Pa/cm. The 
procedures to obtain fiber optic hydrophone sensitivities 
(single and dual coil designs) are discussed in Chapter IV. 
A thorough presentation on fiber optic acoustic sensitivity 
is given by Davis, et al . , CRef. 63. 

C. FIBER OPTIC ACOUSTIC SENSOR 

1. Acousticall y I nduced P hase Modulation 

Laser light conducted by optical fibers submerged in 
a liquid medium may be modulated (intensity or phase) by 
acoustic pressures. This phenomenon is exploited to make an 
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acoustooptic sensor o-f underwater sound by demodulating the 
light using a speci-fic detection scheme. Carome and 
Satyshur CRe-f. 81 present in detail general -fiber mode 
theory and the theory o-F phase and intensity modulation 
induced in optical -fibers. Phase modulation occurs when the 
instantaneous phase angle o-f the lightwave in the optical 
-fiber is varied proportional ly in accordance with the 
instantaneous value o-f a modulating signal, e.g. a 
sinusoidally varying acoustic -field. 

Two major detection schemes -for inter-f erometric type 
sensors discussed in the re-ferences are homodyning and 
heterodyning CRe-f. 5 and 61. Homodyne detection in optical 
inter-f erometry implies that the re-ference beam has the same 
optical -frequency as the sensing beam. On the other hand, 
heterodyne detection is based on the use or mining o-f two or 
more -frequencies, i.e., the -frequency o-f the re-ference and 
sensing beams di-f-fer by Acj. Inter-f erometry is used in 
either case to recover the acoustic signal -frequency. 
Homodyne detection was used throughout this experiment so 
only it will be considered in detail. 

An external pressure -field applied to a -fiber will 
cause various changes in the -fibers physical 
characteristics, such as changes o-f the core radius and 
length, and changes in the values o-f the optical indices o-f 
re-fraction o-f the core and cladding CRe-f. 3 and 41. For the 
case o-f pressure induced changes o-f index and o-f length, the 



30 



phasB shift A, produced by a change of pressure might be 

given asi 

A » k pj^C<0n/3p) + n/J^< P) 3 

o 

(2. 17) 

where n is the optical index of refraction of the core, k 

o 

is the propagation constant of light in the fiber, P is the 
acoustic pressure, and P is the length of fiber subjected to 
the pressure. The partial derivative <dn/0P) is evaluated 
at constant entropy. Whatever the form of the bracketed 
term of equation (2.17), the phase shift amplitude A enters 
the analysis in identical ways so that it is possible to 
proceed without being concerned with its specific form. 

Using a single frequency coherent laser source, the 
time variation of the electric field vector of the 
lightwave, propagating in a fiber that has been subjected to 
an acoustic pressure field, may be expressed ass 

E(t) » E BxpCJ (CJ t -t- Asin( (J t) 3 
o o s 

(2. 18) 

where CJ i* the angular frequency of the coherent laser 
o 

source, CJ * '^1^* sngular frequency of the sound field, and A 
s 

is the phase shift amplitude as described above. 

Phase modulation cannot be directly detected by 
photodetectors due to the fact that they cannot follow the 
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inmtantancous values o-f such high -frequancy variations 

14 

( approx imatsly 10 Hz). Thus« in order to accomplish phase 
detectiony an interf erometric technique must be used to 
convert the phase modulation to intensity (amplitude) 
modulation prior to detection. A monochromatic laser beam 
is split, with part o-f the beam being transmitted into a 

re-ference -fiber. The other beam is sent into a sensing 

-fiber exposed to an acoustic -field. Pressures exerted on 
the sensor -fiber as a result o-f the sound energy density 
produce optical path length changes in the -fiber which, in 
turn, cause phase modulations in the lightwave. A-fter 
passing through the -fibers, the two beams are then 

recombined. Thus the total electric -field is given byi 

~E -1e <t) +“e (t) 

T 1 2 

(2. 19) 



where 



and 



E (t) - E expCJ(CJ t Asin(CJ t) ] 
1 1 o s 



( 2 . 20 ) 



E (t) * E expCJ ( CJ t + cb ) '2 
2 2 o 



( 2 . 21 ) 
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wh®r« 0 



is an arbitrary phase di-F-Ference which will vary 



with temperature, the optical path di-F-Ference, and other 

-Factors. E (t) is the electric -Field vector o-F the sensing 
1 

arm and E (t) is the electric -Field vector of the reference 



arm. The intensity (I) of the recombined beams is 

proportional to the magnitude of the square of E . 



Therefore, after mathematical manipulation as shown by 

Carome and Satyshur [Ref. 8] and neglecting terms that vary 

at angular frequency (jJ *nd 2 6J t since as stated earlier 

o o 

they are undetectable by a photodetector, the following 
equation resultsi 



( 2 . 22 ) 

— k 

where E • E makes explicit the fact that the polarization 
1 2 

directions may not be the same. The first three terms are 
considered the DC terms and all other terms of the series 
are the AC terms. Thus the resulting intensity function 
consists of a series of harmonics of the acoustic 

frequencies where the amplitude of each successive harmonic 
is a function of the acoustic pressure, and varies as the 
Bessel function of corresponding order. 



2 



T 
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If both of the beams are phase modulated by the 

acoustic pressure as in the case of the pressure gradient 

hydrophone, with phase shifts of A and A , respectively, 

1 2 

the final expression for the intensity is identical to 

equation (2.22), except that for A, the quantity 

^A s A - A must be substituted. 

1 2 

The recombined lightwave with its varying intensity 
is then projected onto the surface of a photodetector. 
There results in the photodetector current a signal which is 
related to the acoustically induced phase shift. The 

photodetector current has components of the following form: 

OO 

i(t) = i cos^* fJ (k X) + 2^J (k x)cosC2n((4jt) ]> 

0 0 o / 7 — / 2n o s 

OO 

-i sin(/)* f2 y J (k x)sinC <2n+l) ((J t) 1> , 

0 ^ /)4o2nH-l o 

(2.23) 

where J is the Bessel function of order n CRef. 83,^ is 
n 

the static phase difference, which will vary with 

temperature, path length difference, and other factors, 

k <=2 TT/ \ , where k is the optical propagation constant in 
o o o 

the fiber, \ is the optical wavelength in the fiber, and x 
o 

is the amplitude of the optical path-length change at CJ • 

s 

The J term in equation (2.23) is the DC component 
O 

associated with the acoustically induced signal with the 
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remaining terms being the AC component. Only the AC terms 
of equation (2.23) will be considered in succeeding 
sections. 



2. Fe e dba ck Stabi 1 iz ation 

The purpose of homodyne detection in the fiber optic 

interf erometer discussed above is to transform the optical 

output of the interf erometer into an electrical signal 

proportional to the amplitude of the relative phase shift 

produced by a pressure oscillation at frequency (jJ . 

s 

Homodyne signal processing suffers from phase drift due to 
small temperature and pressure variations. Any random 
change in length in one arm of the fiber interf erometer with 
respect to the other arm appears in equation (2.23) as a 
change in the angle Cf) • As changes, the amplitudes of the 
AC components of equation (2.23) varies due to the factors 
cos and sin(^. This unwanted change in the amplitude of 
the detected signal is referred to as signal fading. 
Consider the case in which the outputs of the two arms of 
the interf erometer are initially in phase, (^ is equal to 
zero, and when combined at the output coupler they interfere 
constructively. If now there is a small acoustically 
induced phase modulation e.g. at peak to peak amplitude 0.1 
radian I referring to Figure 2.S, the photodetector current 
will oscillate about its maximum value at twice the acoustic 
frequency. If the length of one arm increases due to a 
thermal change, the average output intensity of the 
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Figure 2.5 Photodetector Output Current as 

Function of Phase Difference between 
Interferometric Arms 
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interf erometer decreases^ but the acoustically induced 
oscillations will increase in amplitude and approach a pure 
sinusoidal o-f -frequency equal to the acoustic frequency. In 
a practical hydrophone it is essential to maintain the 
interf erometcr continuously at its maximum sensitivity. The 
sensitivity to phase changes varies as the slope of the 
curve of Figure 2.5. The maximum sensitivity occurs for odd 
multiples of 7T/2 while essentially zero sensitivity occurs 
for even multiples of TT/2. Therefore to maintain 
sensitivity a 7T/2 phase bias is required. This 7T/2 bias 
condition is known as quadrature. 

To place the interf erometer in quadrature condition 
for maximum sensitivity and in addition provide a means to 
compensate for slow drifts, a phase-locked compensator 
circuit is used. The circuitry required to provide and 
insure quadrature operation in the presence of low frequency 
drift is shown schematically in Figure 2,6 CRef. 91. The 
outputs of the two photodiodes are combined in a 
differential amplifier that provides common-mode rejection 
as well as amplification. This is followed by two stages of 
integration that provide additional amplification. These 
two integrator — amplifiers pass all signals from DC up to the 
highest frequency of interest. The output of the two stage 
integrator — amplifier is applied to a phase shifter located 
in the reference arm of the interf erometer . The phase 
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Figure 2.6 
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Phase-Locked Loop Homodyne Detection 
Circuit 
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shifter consists of a lead zirconate-lead titanate 
piezoelectric (PZT) cylinder around which the fiber in the 
reference arm is rather tightly wound. The output of the 
integrator — amplifier is Just equal to the voltage required 
to compensate for phase changes associated with the low 
frequency drift and that induced by the acoustic pressure 
signal of interest. The effect then is to produce a phase 
shift in the reference arm almost equal to that in the 
sensing arm, causing the interf erometer to remain very close 
to perfect balance, i.e. , to phase-lock the system CRef. 
103. If the phase were exactly locked there would be no 
output signal from the interf erometer . However, there must 
be an error signal or resulting phase modulation at the 
photodetectors in order to have a feedback signal and this 
is effectively 1/G times the modulation that would exist 
without the feedback, where G is the overall gain of the 
feedback circuit. If the system is initially at a bias 
(operating or quiescent) point away from quadrature, the 
effective gain will not be at its maximum value and the 
system will tend to drift toward its maximum gain point and 
therefore toward quadrature. The output signal of the 
feedback system represents the correction signal required to 
keep the phase difference between the interf erometer at 7T/2 
radians. It contains drift noise and signal information. 

The signal out of the compensating circuit is also 
fed through a high-pass filter that has its low frequency 
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limit sBt at the lowest frequency of interest but above the 
frequency of thermal drift etc. Therefore, the resulting 
output is a band of frequencies corresponding to the 
frequency range of interest. This constitutes the output of 
the interf erometric sensor. 

Operational amplifiers <OPAMPS) are used in the 
feedback circuit and complementary metal oxide semiconductor 
(CMOS) components in the reset circuit. The levels of 
voltage that can be applied by these circuits to the phase 
shifter are the order of ± 10 volts. On the other hand, the 
range that the piezoelectric phase shifter can accommodate 
is hundreds of volts. Furthermore, in many cases the 
amplitude of the phase drift resulting from temperature or 
pressure changes is much larger than the phase shift that 
would be generated by applying ± 10 volts to the phase 
shifter. Thus, the circuitry keeps track of how much 
voltage has been applied to the phase shifter and if the 
limit of the circuit begins to be reached it automatically 
and rapidly resets the circuit back to its mid range from 
which point it restarts the compensation process. This is 
the purpose of the reset circuit indicated in Figure 2.6. 
The phase change associated with a large amplitude slow 
drift is compensated by a number of saw toothed— like small 
amplitude phase changes. Care must be taken in the 
alignment of the compensation circuit to minimize noise, 
introduced during the reset process. 
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III. EXPERIMENTAL APPARATUS 



A. INTERFEROMETER DESISN AND CONSTRUCTION 

Figura 3.1 shows a schematic diagram o-f the experimental 



interf erometric 


system 


utilized 


in this 


study. 


The 


interf erometer 


is modular in construction. 


As each 


module 


was completed 


proper 


operation 


was ver i f i ed . A 


brief 


description of the 


overall 


construction of 


the 


interf erometer 


is first 


presented) 


this is 


followed 


by a 



detailed description o-f each individual part. 

The system is comprised o-f an all —fiber Mach— Zehnder 

inter-f erometer . A Laser Diode Labs solid state laser was 
1 

pigtailed to high numerical aperture single-mode (830 nm) 
Corning -fiber. The laser, with the pigtailed -fiber, was 
spliced (see section on -fiber preparation and splicing) to 
approximately hal-f a meter o-f ITT single-mode fiber. The 
free end of the ITT fiber was then spliced to an ITT 2X2 
all-fiber 3dB coupler. Initial attempts were made to splice 
the pigtailed Corning fiber directly to the coupler fiber 
but with limited success. Differences between melting 
temperatures 'of the Corning fiber and coupler fiber resulted 



1 Pig tailing of fiber to laser accomplished by laser 
manuf acturer 
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I. Laser diode 5. Sensing coil 

2.3dB Coupler 6. Phofodefecfor 

3, Polarization controller 

4. PZT Phase shifter X-Indicates splice 



Figur* 3.1 Schsmatic o-f Exparimantal 
Intar-faromatric Systam. 
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in low strength, high loss splices. For this reason the ITT 
T-1601 fiber had to be used as an intermediate link between 
the Corning fiber and ITT coupler fiber. 

The couplers themselves are all fiber counterparts of 
bulk optical beam splitters and equally split the light at 
the input end of the interf erometer . Half the light is 
coupled to one arm and the remaining light into the other 
arm. 

Two single~mode optical fibers with approximately equal 
length (estimated path length difference of approximately 2 
cm) form two arms of the interf erometer . Both were 
constructed identically; therefore only the construction of 
one will be discussed. Eight meters, from a ten meter 
length of fiber, was wrapped in a single layer onto a 
piezoelectric (PZT) cylinder. At each end of the wrapping 
approximately a meter of lead was left. One lead was then 
spliced to fifteen meters of fiber of which approximately 
ten meters was coiled into a 4 cm diameter coil. The other 
lead coming from the fiber wrapped around the PZT is fed 
through mechanical polarization controllers CRef. 111. One 
end is then spliced to one of the output leads of the input 
coupler and the other end is spliced to one of the input 
leads of the output coupler to complete construction of an 
arm of the interf erometer . The output coupler recombines 
the light from the two arms. Basically the interf erometer 
is used to convert phase modulation to intensity modulation 
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which occurs when light ■Froin the two arms is combined. The 



two outputs of 


the 


i nterf erometer 


are 


detected 


by two 


photodiodes. 












Figure 3.2 


is 


a photograph 


of 


the completed 


i n t er f er omet er . 


The 


design of the 


system allows 


for a 



single coil to be placed in water to act as a pressure 
hydrophone or for both coils to be placed in water for 
operation as a pressure gradient hydrophone. 

B. INTERFEROMETER COMPONENTS 

1. Laser Sour ce 

The optical source used in this experimental 
interf erometric system was a single frequency infrared laser 
diode optimized to emit coherent light at a wavelength of 
830 nanometers. Two solid state laser diodes from Laser 
Diode Labs were used. The first was a CW-234 laser, 
pigtailed by the manufacturer with a short section of high 
numerical aperture single mode (830 nm> Corning fiber. The 
second laser, a CWX-233 laser also, pigtailed with Corning 
fiber, replaced the first laser after an inadvertent 
increase in current destroyed the CW-234 laser. Power for 
the laser was provided by a 12 volt 18 amp-hr Gel cell 
battery through a voltage regulating circuit. This circuit 
maintained constant current to the laser diode. Figure 3.3 
is a schematic of the power supply. Before using either 
laser in the interf erometer the end of the Corning fiber was 
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Figure 3.2 



Photograph of Experimental 
Interferometric System 
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Figure 3.3 Schematic of Laser Diode Power Supply 



46 



cleaved and set up to enable a Coherent 212 optical power 
meter to sense the laser diode power output. As current to 
the laser diode was increased -from 0 to approximately 53 
milliamps -for the CW-234 laser, and 0 to 39 milliamps -for 
the CWX-233 laser, the power o-f light emitted -from the laser 
and coupled into the -fiber was read on the power meter. 

Figures 3.4 and 3.5 are graphs o-f laser power vs source 

current -for the CW-234 and CWX-233 lasers, respectively. 
Extrapolating the linear region back to zero power yields a 
threshold current of 48.5 milliamps for the CW-234 laser and 
31.3 milliamps for the CWX-233 laser. The current to the 
laser diode was maintained at S3 milliamps for the CW-234 

laser and 35 milliamps for the CWX-233 laser during 

experimental testing and data acquisition. 

2* §E®£lfi£§tl9Q 

The optical fiber used was ITT single-mode fiber, 
2 

type T-1601 optimized at 0.83 microns. Specifications of 
fiber are as follows: 



Fiber I dent: 

Preform No. : 
Numerical Aperture: 
Core Diameter: 



81092U-17Bic 
EMT-21972B 
0. 12 

4.57 microns 



2 Optical fiber obtained from International Telephone 
and Telegraph, Electro Optics Products Division, Roanoake, 
Va 24019 
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225 -L LASER DIODE LABS CW-234 SOLID STATE LASER 
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LASER DIODE DRIVE CURRENT (MILL I -AMPS) 
Figure 3.4 CW-234 Laser Diode Power Curve 
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Outer Cladding Diameter i 75 microns 



Primary Sheath! 


GE 615 Silicone 


Secondary Sheath: 


polyester Hytrel 


Total Diameter: 


16 mil 


Attenuation: 


2.07 dB/Km at 0.83 microns 



Figure 3.6 is a sketch, provided 


by the manufacturer of the 



end of the fiber, with typical dimensions, approximate index 
of refraction profile, and typical spectral attenuations. 

3. 3 dB Couple r 

The 2 X 2 single-mode couplers were also purchased 
from ITT. The following specifications are provided! 



Serial Nos.: 


JM-SM-108 & JM-SM-109 


Fiber No.: 


EMC 41556C/830427-401a 


Fabrication Date: 


6/3/83 


Excess Loss: 


0.1 dB 


Uniformity! 


0.1 dB 


Operating Wavelength: 


0.83 microns 



The purpose of a 3 dB coupler is to split the light equally 
into the two arms of the interf erometer or to recombine the 
light causing the light from the two arms to interfere in 
the output fibers and on the face of a photodiode. 
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AnENUATION (dB/Kin) 



DIMENSIONS SHOWN ARE NOMINAL VALUES 




R 



, COPE ^ SILICA 
n CLACOING 



I 



INDEX Of REFRACTION PflOFiLE 




TYPICAL SPECTRAL ATTENUATION - SINGLE MODE OPTICAL FIBER 

Figure 3.6 ITT Single Mode Fiber T-1601 
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4. PiezoBl_ectric P ha sB ShlftBr 

Tha phasB shi-ftar consists of a lead zirconate-lead 
3 

titanate (PZT) cylinder around which the fiber is tightly 
wrapped. The cylinders used were type Channel 5500. Each 
cylinder is 3.8 cm (1.5 in.) long by 3.8 cm in outer 
diameter with a wall thickness of 0.32 cm (0.225 in.). Two 
of these cylinders were incorporated into the 
interf erometer , one in each arm. The PZT cylinders have a 
relatively small dynamic range but by wrapping many turns of 
fiber on the PZT it was possible to increase the magnitude 
of the phase shift. Chapter II, Section C presents the 
theory of operation of the PZT as a phase shifter. 

5. Pol ar izati on Cont ro l ler 

A polarization controller, as described by Lefevre 
CRef. 113, was used in this experiment. A very simple 
configuration is employed to control the polarization of 
single mode fiber. This device is equivalent to fractional 
wave plates of classical optics. The controller uses the 
stress birefringence induced by bending the fiber. 

Photodete ct or 

The photodetector used to detect the optical output 

4 

of the fiber interf erometer is a Clairex type CLD-41 



3 Piezoelectric cylinders obtained from Channel 
Industries, Inc., Santa Barbara, California 93105 

4 photodiodes obtained from Clairex Electronics, Mount 
Vernon, New York, New York 10550 
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photodiode. They are all silicon PN planar diodes with high 



linearitVf low dark current, and -fast response. The 
•following are the electrical characteristics as given by the 
manu-f acture: 



Active Area: 

Short Circuit Currenti 
Open Circuit Voltage: 
Dark Current: 

Junction Capacitance: 
Rise or Fall Time: 



1 . 3x 1 . 3 mm 

min 6 to max 12 microamps 
.40 volts typical 
1 microamp 
200 pi co-farads 
5 microseconds 



Temperature Coe-f -f icient : +. 2%/degree C typical 
Peak Spectral Response: 0.91 microns 

Figure 3.7 is a plot o-f photodiode DC voltage vs system 
laser optical power as measured in the laboratory using a 
Hewlett-Packard HP-3438A digital multimeter to measure the 
output voltage. 



C. FIBER PREPARATION AND SPLICING 



Be-fore splicing between -fibers could be accomplished, 
each end o-f the -fibers had to be prepared by stripping and 
cleaving. Fiber stripping was accomplished in two steps. 
First, the outer Hytrel plastic sheath was removed by use o-f 
a razor blade. The fiber was placed on a piece of black 
construction paper and then a razor placed at a very small 
angle to it. With a slight pressure applied to the razor, 
the fiber was drawn out from underneath the razor. This 
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37± CLAIREX ELECTRONICS CLD-41 PHOTODIODE 
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LOG RELATIVE OPTICAL OUTPUT POWER-MICROWATTS 
Figure 3.7 Clairex CLD-41 Relative Irradiation Curve 



will carve the Hytrel jacket neatly -from the fiber. By 
rotating the fiber and repeating this process several times, 
most of the outer Hytrel and inner silicone sheathing 
material will be removed. 

The second step, that of removing any residue of inner 
or outer sheath that was not removed by the razor blade, is 
best done chemically with a solution of chromic acid. This 
is done by immersing the section of fiber that has been 
partially stripped by the razor blade in the acid for 
approximately two minutes. Any Jacket material that remains 
will turn to a Jelly-like consistency which is easily 
removed by wiping with a soft tissue. After rinsing in 
distilled water, the fiber is pulled between a tissue 
saturated with methyl alcohol. If the fiber is completely 
clean of the Jacket a squeaking sound will be heard. At 
this point the fiber is bared to the silica cladding and 
must be handled with caution. An unscratched fiber has very 
high tensile strength, but will break easily if scratched 
and bent. 

In order to insert and extract light efficiently to 
obtain a low loss splice, the ends must be flat and 
perpendicular to the longitudinal axis of the fiber. The 
process utilized to achieve this was cleaving. To cleave 
the fiber, the ends are first stripped back about two 
centimeters longer than needed. After preparation for 
cleaving, the fiber was held with slight tension and 
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carefully touched to the cleaving blade at a point 

5 

approximately 2 cm from the end of the fiber. The cleaved 
area was then placed on a smooth, round mandrel that has a 
radius of approximately four centimeters. The fiber was 
held with thumb of one hand. Maintaining modest tension, 
the fiber was pulled around mandrel with other hand. When 
scored area reaches mandrel, a clean break will usually 
result. Experience showed that if the fiber broke while 
being scored it would have a bad end face. 

To check the cleave, the fiber is inspected under a 
microscope. Angular breaks other than 90 degrees will not 
make a good low loss splice. In addition, chips and barbs 
at the cleaved ends will not allow satisfactory splices. 
Figure 3.8 shows a good cleave along with improper cleaves. 

Splicing of fiber was accomplished using an electrical 
arc welder. The type of welder used was a Power Technology, 
Inc. portable fiber splicer model PFS-240. Basic splicer 
operations of prefusing and fusing for monomode fiber were 
done in accordance with the operating manual. Thumbwheel 
settings which controlled splicer arc amperage and the time 
variation of the arc current usually used are as follows: 
Ramp Time Arc Time Arc Current 
0.1 sec 0.6 sec 16 mA 



5 Sapphire cleaver obtained from General Fiber Optics, 
Caldwell, New Jersey 07006 
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Figure 3.8 Representation of Proper and Improper Cleaves 



Final alignment of cores -for maximum light transmission 
was used to produce a low-loss splice. The 
micro-positioners on the splicer ware adjusted to bring the 
prepared ends o-f the -Fiber into close alignment. Light -From 
the laser was injected into the -Free and o-F one o-F the 
-Fibers, trans-Ferred across the proposed splice to the second 
-Fiber, and the output -From the latter was projected onto the 
-face O-F the power meter's photodiode. The micropositioners 
were adjusted until the power meter indicated that the 
transmission across the splice was maximum. At that point 
the -Fiber ends were brought together, heated by the arc 
welder and -Fused. As the Joint cools, sur-Face tension will 
pull the -Fibers into -Final alignment. Reports in the 
literature indicate that losses obtained by -Fusion splicing 
may be as low as 0.1 's dB. However, in this study many 
hours O-F practice were required to obtain good (i.e. less 
than 3dB loss) splices. Nina splices are required in the 
con-Figuration o-F the experimental system described above. 
The average loss -From each splice was approximately IdB and 
the total loss o-F the constructed inter-Ferometer was 
approximately 8.4dB. 

Prior to projecting light -From the end o-f the -Fiber onto 
the sur-Face of a photodiode, the end of the fiber was 
prepared in the same manner as for splicing with the 
exception that the length of fiber that had its sheath 
removed was painted black. Painting was used to accomplish 
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mode stripping. Cladding light was stripped prior to 
reaching the end of the -Fiber, thus only core light was 
detected by the photodiode. In splicing each additional 
-fiber element into the system it was necessary to mods strip 
its output end Just be-fore the energy meter to ensure that 
only core light is detected, i.e. to absorb unwanted light 
that may reach the output end via the cladding. 

D. INSTRUMENTATION AND DATA ACQUISITION 

Basic block diagrams o-f the experimental systems used to 
obtain data on intsr-ferometer per-formance will be shown in 
Chapter IV as the experimental procedures are discussed. 
Computer data acquisition was used -for portions o-F the data 
taking. The Hewlett*-Packard HP-85F computer coordinates the 
peripherals and records and displays the data. The printer 
and plotter are non-essential to the system and were used to 
provide larger -format graphics and printouts than are 
available with the built-in thermal printer provided with 
the HP-8SF. The -following is a brie-f description of each 
individual piece of equipment. 

1 • bPz§5E Comput er 

The HP-85F is an eight bit microprocessor that 
utilizes BASIC computer language. The computer has as 
standard 16K bytes of read/write memory. 16K bytes of 
additional memory were added to the system giving a total of 
32K bytes. The HP-8SF has a 127 millimeter diagonal black 
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and white electromagnetic CRT. A 32 character per line 
thermal printer/plotter is part of the unit. Programs or 
data may be stored on and read from magnetic tape 
cartridges. To interface with peripheral equipment, an I/O 
ROM and an interface card were added to provide HP-IB (IEEE 
standard 488-1975) instrumentation capabilities. 

2. Arbitrary W av eform Generator HP 33^40 

The Hewlett-Packard Model 3314A arbitrary waveform 
generator can produce three different pre-programmed 
waveforms! sine, square, and triangular. This instrument 
has a frequency range of 1 microHertz to 20 megaHertz. 
Output amplitude is 0 millivolts to 10 volts peak to peak 
into a fifty ohm load. This model is fully programable 
through a HP-IB connection. For this experiment sine waves 
of various frequencies and amplitudes were usually 
generated. During alignment of the locking circuit, 
however, a small amplitude 100 Hertz square wave was used. 

3- Sgectrums A na ly zer HP 3S82A 

The Hewlett-Packard Model 3582A is a dual channel 
spectrum analyzer. This instrument has a frequency range of 
0.02 Hertz to 25,600 Hertz. The analyzer has a 11.9 by 9.6 
cm CRT that can display two simultaneous information traces, 
plus four lines of alphanumeric data giving measurement 
configuration and results. Frequency spans from 1 Hertz to 
25,000 Hertz full scale allow flexibility in selecting the 
portion of the spectrum to be analyzed. Spans from 5 Hertz 
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to 25,000 Hertz can be positioned anywhere within the 
-Frequency range o-f the instrument to provide excellent 
-Frequency resolution. The instrument's “-Front-end" 
sensitivity ranges can measure and analyze -From 1 micro-volt 
to 31.6 volts and has a dynamic range o-F 70 dB. 

4. ?215 

The Tektronix Model 2215 oscilloscope is a dual 
channel, 60 MHz instrument. Its vertical system provides 
calibrated deflection factors from 2 millivolts to 10 volts 
per division. Trigger circuits enable stable triggering 
over the full bandwidth of the vertical system. The 
horizontal system provides calibrated sweep speeds from 0.5 
seconds per division to 50 nanoseconds per division. 

5l9l^§l Multimet er HP 3438A 

This multimeter was utilized to obtain rms voltage 
readings of an LC-10 hydrophone output and transmit this via 
the HP-IB to the HP-85 computer. 

6. Standard Hydrophone LC-1,0 

6 

An LC-10 hydrophone , Serial No. 2341 was used as 
the standard hydrophone for sensitivity calculations of the 
single coils and of the dual (gradient) coil setup. Average 
free field voltage sensitivity for this hydrophone is 
specified by the manufacture to be -209.2 dB re 1 Volt//iPa. 



6 Hydrophone manufacture is Celesco Ind. , Canoga Park, 
California, 91304 
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This hydrophons and an USRD type F27 reversible transducer 
were used in another experiment CRef. 121 to obtain a 
reciprocity calibration of an USRD type F33 transducer. The 
sensitivity for the LC-IO hydrophone obtained from the 
experiment was within 2 dB of manuf acturer 's specifications 
between 5 an 30 KHz. 

7. Acous ti c Calib ra tor 

An acoustic system was constructed for calibrating 
hydrophones of various types by the comparison method. 
Figure 3.9 shows the calibrator in the lab, and Figure 3.10 
shows the details of construction. The calibration section 
consists of an open-ended column of water excited by a 
moving-coil driver. A USRD type Jll transducer is used as 
the driver. The calibration chamber consists of a 6 inch 
inner diameter, 1/4 inch thick PVC 1120 tube (A8TM 1785), 
mounted on a rim around the face of the driver. 
Compensation for the hydrostatic pressure of the water on 
the face of the driven piston is provided by the airtight 
enclosure mounted under the driver assembly. A manometer is 
used to measure the equalizing pressure. A valve is opened 
and air is pumped into the equalizing chamber until the 
internal air pressure is equal to the water pressure on the 
upper face of the driver. This system is similar in design 
to the USRD type 840 shipboard calibrator used by the U.S. 
Navy CRef. 133, based on article by Schloss and Strasberg 
CRef. 143. 
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Figure 3.9 



Photograph of Acoustic Calibrator 
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Figure 3.10 Sectional View o-f Acoustic Calibrator. 
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The compliance o-f the PVC tube reduced the speed of 
sound of the water in the calibrator to approximately 1/5 of 
the normal speed of sound in water. This was useful in the 
experiments since it reduced the wavelength X by a factor 
of 5 at any given frequency and resulted in a larger 
gradient. Figures 3.11, 3.12, and 3.13 show the standing 
wave acoustic field at three resonant frequencies of the 
calibrator. The figures are graphical plots of LC-10 
hydrophone voltage as a function of position in the 
calibrator. By measuring the distance of a half wavelength 
at the three resonant frequencies the speed of sound in the 
calibrator can be calculated. From the data shown in Table 
I an average value of 358.3 meters per second with a 
standard deviation of 8.6 m/sec for the speed of sound in 
the calibrator is calculated. The data was gathered over a 
period of one month. 



TABLE 



Sound Speed in 



Frequency (Hz) 


Half 

Wavelength 


390 


44.5 


390 


44 


391 


45.5 


595 


30 


593 


31 


593 


30.5 


593 


30.5 


593 


30.5 


752 


24.5 



Cal ibrator 

Speed of 

(cm) Sound (m/sec) 

347. 1 
343.2 
355.8 
357.0 

367.7 

361.7 

361.7 

361.7 
368.5 



average >>> 358.3 -»■ 8.6 
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DISTANCE FROM TOP OF WATER COLUMN-CENTIMETERS 

Figure 3.11 LC-10 Voltage in Calibrator at 391 Hz 



STANDING WAVE PATTERN IN CALIBRATOR AT 593 HERTZ WITH 2 VOLTS DRIVE 
DEPTH OF WATER COLUMN 91.5 CENTIMETERS 
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Figure 3.12 LC-10 Voltage in Calibrator at 593 Hz. 
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DISTANCE FROM TOP OF WATER COLUMN-CENTIMETERS 



5.4- STANDING WAVE PATTERN IN CALIBRATOR AT 752 HERTZ WITH 2 VOLTS DRIVE 

DEPTH OF WATER COLUMN 91.5 CENTIMETERS 
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DISTANCE FROM TOP OF WATER COLUMN-CENTIMETERS 

Figure 3.13 LC-10 Voltage in Calibrator at 753 Hz 



8. Lo ck ing Ci_rcui_t 

A phasa lock loop locking circuit similar to the 
circuit presented by Fritsch and Adamovsky CRef. 91 and 
discussed in Chapter II was used to verify operation of the 
interf erometer as a single and dual (gradient) hydrophone 
system. Data will be presented in the next chapter. 
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IV. EXPERIMENTAL PROCEDURE AND RESULTS 



A. INTERFEROMETER CHARACTERISTICS 

A-fter construction of the interferometer was completed « 
testing was conducted to verify proper operation of the 
interf erometer as a sensor. The initial tests consisted of 
detailed measurements taken to determine the amplitude of 
the optical phase shift as a function of the drive voltage 
applied to the piezoelectric cylinder in one of the arms of 
the interf erometer . Figure 4.1 shows a block diagram of the 
instrumentation required for operating the system. A 
variable amplitude sine wave of frequency f from the 
function generator was applied to the PZT cylinder. A 
Tektronix 2215 oscilloscope and a HP 3S82A spectrum analyzer 
were used to detect the output of one of the photodiodes. 

The behavior of a Mach-Zehnder fiber interf erometer was 
treated in detail in Chapter II. Recall that the AC portion 
of the photodiode current is given by the following 
equation: 



i (t) > i cos 
0 



r ^ 

0 *[2 g, 



J (k X )cosC2n ( Cl) t ) 1 
2n o s . 



Ln0 



-i sinffi* 2 / J 



(k x)sinC (2n*t-l) ( CJ t) 3 



n=0 2n+l o 



(4.1) 



where J is the Bessel function of order n. 
n 
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Figure 4.1 Block Diagram of Instrumentation Setup 
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From equation (4.1) it is clear that as the phase shift 

k X is increased the photodiode output becomes non-linear 
0 

and harmonics of the drive signal begin to appear. This is 
illustrated in Figure 4.2 where oscillograms of the 
photodiode output signal are shown for successively higher 
drive signals applied to the piezoelectric cylinder. In 
each oscillogram the drive signal is shown in the upper 
trace and the corresponding photodiode output signal is 
shown in the lower one. 

The magnitude of the optical phase shift as a function 

of drive voltage may be determined quantitatively by 

examining the frequency spectrum of the photodiode output 

signal. From equation (4.1) it may be seen that the 

th 

amplitudes of the various n order harmonics should have 

nulls for phase shifts k x equal to the zeroes of the 

0 



th 

corresponding n order Bessel functions. This behavior is 
illustrated in Figure 4.3 where the photodiode voltage is 
plotted for the fundamental and first two harmonics vs drive 
voltage applied to the piezoelectric cylinder. The data 
were obtained with the computer data acquisition equipment 
discussed in chapter III. A computer program for the HP-85F 
computer, listed in Appendix A, was written in BASIC 
computer language to control the peripheral equipment. A 
fundamental frequency, lower and upper limits of amplitude, 
and an incremental value are selected. The function 
generator and spectrum analyzer are preset and configured to 
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/ 
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Figure 4.2 



Oscillograms of PZT Drive Signal and 
Photodiode Output Signal 




a 



the prescribed settings. At the lower voltage amplitude 
limit the spectrum analyzer seeks -for the maximum amplitude 
of the fundamental and first two harmonics attained during a 
selected period of time. The amplitude of these signals 



time selected on the spectrum analyzer, of the order of 16 
seconds, was sufficient that the actual maximum values were 
attained and recorded. The HP-85F computer directs the 
spectrum analyzer to output the maximum amplitude values of 
the fundamental and harmonics. These values are stored by 
the computer. The computer then directs the function 
generator to increase the drive voltage by the specified 
incremental value. The spectrum analyzer again seeks for 
the maximum amplitudes and the process is repeated. This is 
done for each incremental value until the upper limit is 
reached. The data stored in the computer memory is sent to 
the printer and also transferred to magnetic tape for 
permanent storage to enable analysis at a later time. A 
graph of the data is plotted on the computer printer. 
Sample of data obtained is seen in Appendix B. From this 
data plots of the type shown in Figure 4.3 could be 
generated. 

The computer program was run initially with the CW-234 
laser source. Three runs were made, two with the fiber 
wrapped PZT cylinder of one arm acting as the modulator and 
one with the other fiber wrapped PZT of the other arm. 



vary in time as 




varies randomly. The sample 
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A-fter replacement o-f the source with the CWX-233 laser was 

accomplished two more runs were made, one on each PZT 

cylinder. From the data, the amplitude o-f the drive signal 

that causes the -fundamental to be at its first null was 

determined to be 1 . 65V (peak) +0. 04V (peak) . Thus, referring to 

equation (4.1), it may be concluded that the drive condition 

was such that the first order Bessel function J was at a 

1 

zero. The first of the fundamental 's zeroes occurs at a 

phase shift of k x a 3.83 radians. Therefore for the 8 

O 

meters of fiber wrapped on the PZT cylinder a value of 0.29 
rad/volt (peak) /meter results. By observing the drive 
voltages required to produce the first zero of the 
fundamental as the frequency was varied manually it was 
possible to determine that the phase shift per volt remained 
constant, within the precision of the measurements, for 
frequencies from 80 hertz to 1500 hertz. 

Agreement of experimental data presented in Appendix B 
and illustrated in Figure 4.3 to Bessel function theory is 
shown in Table II. The ratio of experimental maxima and 
minima of the same harmonic are compared to the theoretical 
ratio of maxima and zeros of Bessel functions of the same 
order. Theoretical data as tabulated by Abramowitz and 
Stegun CRef. 153 was used. 

The last test conducted was to verify operation of the 
system with the locking circuit. After aligning the 
locking, a small amplitude sine wave was applied to one of 
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Comparison of Experimental Data with Bessel Theory 
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Th» 



the PZT cylinders simulating an acoustic signal, 
locking circuit produced a signal, similar to this input 
signal, that was fed back to the PZT cylinder in the other 
arm to act as a compensation signal, thus maintaining the 
two arms in quadrature. This output signal from the locking 
circuit was also sent to the oscilloscope to show that the 
output signal did follow the input signal. 

B. INDIVIDUAL SENSOR SENSITIVITY 

The individual acoustic sensitivity of the sensor coils 
in each arm of the interf erometer was calculated from 
measurements obtained in the acoustic calibrator. A single 
coil was mounted on a stand with an LC-10 hydrophone placed 
in the center of the coil as shown in Figure 4.4. The 
acoustic calibrator was driven at a frequency that generated 
a standing wave acoustic resonance. Chapter III section D7 
discusses the standing waves generated in the calibrator. 
The coil /LC-10 hydrophone setup was lowered to the first 
maximum point of the standing wave. Output of one of the 
photodiodes of the interf erometer was sent to the Tektronix 
2215 oscilloscope and the HP-3582A spectrum analyzer. 
Output of the LC-10 hydrophone was sent to the HP-3438A 
Multimeter set to read AC voltage. A schematic of the 
instrumentation setup is shown in Figure 4.5. The spectrum 
analyzer was set to monitor the fundamental frequency and 
the first two harmonics. Drive voltage was increased to the 
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Figure 4.4 Photograph of LC-10 Hydrophone and 
Single Fiber Coil Setup 
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Figure 4.5 Block Diagram of Calibration Circuit 
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acoustic calibrator until the voltage of the fundamental 
frequency, as seen on the spectrum analyzer, reached its 
first minimum value. The voltage of the LC-10 was 
recorded. If there was more than one pressure maximum the 
sensor arrangement was moved to the next pressure maximum 
and the LC-iO voltage at the fundamental 's first zero point 
was recorded. Data were recorded for all of the pressure 
maxima of the standing wave. Data were obtained both 
manually and by computer control as discussed in the above 
section. 

The first minimum of the fundamental occurs at a phase 

shift of 3.83 radians, therefore the sensitivity of the coil 

M in optical radians per micropascal is calculated as 
C 

follows: 



M = M *3.83rads/V 
C H H 

<4. 2) 



where M is the sensitivity of the LC-10 hydrophone in 
H 

volts/^Pa obtained from manufacture specification sheets. 

V is the LC-10 hydrophone voltage where the first zero of 
H 

the fundamental occurred. Table III is a summarization of 
all single coil sensitivity measurements taken at three 
different calibrator resonances and different depths 
associated with pressure maxima. 
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Table III 



Individual Coil Sensitivity 



Coil # 1 ! : Coil « 2 

! : 



Freqi Depth of 
(Hz)!Hyd from 
1 Air-Water 
S Boundary 

! (cm) 

1 

1 


Nos. 

of 

Runs 


SAvg 
! M 
iirad 


Std ilFreq! Depth of 
dev i!(Hz):Hyd from 
A^ad!! lAir — Water 


: Nos. 
: of 

: Runs 

1 

1 

1 

1 

: 


:Avg ;std 
: M : dev 
:/irad :4Zr ad 


ill Pa 

• -3 
1 10 


^p.i: 

;; 


Boundary 

(cm) 


'./i PaiifjL Pa 
:io‘^ :io'^ 


392 1 


23 


1 


17.07 


- \ S390 1 


! 23 


: 1 


:5.65: - 


385 ! 


60 


1 


: 12.4 


- \ 1 389 1 


1 65 


: 1 


:9.99: - 


385 1 
1 


80 

1 


9 

! 


19. 15 
1 


2. 31 :: 385 i 78 

1 1 1 


: 9 
1 


: 10 . 3 : 2.4 

1 1 


1 

544*1 


( 

16 1 


1 

1 6 


1 

14.95 


1 1 1 
.563: :593 1 


1 16 


1 

: 4 


1 1 

:5. 19 : .314 


593 i 


48 


7 


!8.66 


1.14: :593 1 


1 48 


: 3 


: 1 1 . 5 : . 529 


544 ! 
1 


47 


10 


19.93 

1 


2.97: : - 
1 1 


— 


1 _ 
1 

! 


• « • _ 
1 1 

1 1 


1 

695* t 


12 


9 


1 

19.77 


1 1 

2.41 : :704 


11 


1 

: 13 


1 1 
'• 8 . 45 : 1.36 


695 : 


37 


5 


: 10.4 


. 295 : : 695 


36 


: 18 


: 12.6: 1.57 


695 ! 


65.5 


7 


!25.4 


4.99: :700 


62.5 


: 18 


:2i.6:9. 15 


695 ! 


91.5 1 


1 28 


I18.6S6.8 ::690 


92.5 


: 6 


: 20. 7 : 9. 27 



* Depth o-f water column changed from 93cm to 100cm to 
shift second resonance from 593 to 544 Hz « and to 
shift third resonance from 752 to 695 Hz. 



If all of the measurements taken are grouped together 
neglecting frequency and coil position in the calibrator, 
the blind average of the sensitivities are: 

-3 

Coil #1 M = 13.8 X 10 /irad//iPa 
Cl 

with a standard deviation of: 



-3 



0" = 7.5 X 10 ^rad//iPa 



82 



and 



-3 

Coil #2 M = 13.9 X 10 iL/ rad//i Pa 
C2 

with a standard deviation o-ft 

-3 

(J = 7.7 X 10 ^ rad/^Pa 

As can be seen -from Table III and the blind averages, the 
individual coil sensitivities are matched well within 
experimental error. 

Table III indicates an increase in coil sensitivity as 
depth is increased with exception of very deep depths where 
proximity of the sensors to the source appears to have 
effected the measurements. 

C. GRADIENT SENSOR SENSITIVITY 

A gradient sensor was formed by placing the two fiber 
coils in the calibrator separated by a distance of 10 cm. 
The LC— 10 hydrophone is placed in the center of the coils 
approximately 5 cm from each coil. The LC-10 hydrophone and 
both of the coils are then submerged in water and exposed to 
an acoustic field and an established acoustic pressure 
gradient. If the pressure sensitivities of the two coils 
are closely matched, the pair does not respond to the 
pressure, but to the spatial pressure gradient. Such a 
gradient sensor can be utilized to detect the direction of 



83 



an incident acoustic wave as discussed under conventional 



gradient hydrophones in chapter II. 

The Jll transducer source o-f the acoustic calibrator is 
driven at a frequency to produce a well defined standing 
wave. The LC— 10 hydrophone/ dual coil setup is submerged to 
a location of a pressure minimum, defined as a point where 
the voltage from the LC-10 hydrophone is at a minimum. With 
the LC-10 at this position the two coils are located on 
either side of the pressure minimum and therefore in a 
region of maximum pressure gradient. The external 
instrumentation arrangement is the same as when individual 
coil sensitivities were obtained. The driver amplitude is 
increased until the amplitude of the fundamental frequency 
reaches its first zero point. However, the drive voltage, 
not the LC— 10 voltage was recorded. Since the LC-10 
hydrophone is at a minimum pressure position, its output 
voltage is quite low and, therefore, the pressure in the 
region of the coils could not be directly computed using 
just LC— 10 sensitivity. 

However, by knowing LC-10 voltage (V ) at a pressure 

H 

maximum for a set drive voltage a correlation between LC-10 

voltage and Jll drive voltage could be calculated. Using 

LC-10 sensitivity (M ) , a value of pressure per Jll drive 

H 

volts is obtained. This is shown in the following formula: 

P /Jll volts = CV at P D/CM «JllvoltageD 
o H max H (4.3) 
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Recall front Chapter II, equation (2.6), the pressure at 
a distance x above and below (*i'|-) a pressure minimum 
(assuming x <<X) is: 

P = +P kx 
+ o 

(4.4) 



Replacing P of equation (4.4) with equation 
o 

multiplying by the drive voltage needed to reach 

zero of the fundamental (V ) the pressure at a 

Jll 



is: 



(4.3) and 
the first 
distance x 






C (P /Jllvolts)*V Ikx 
o Jll 



(4.5) 



With the coils equally spaced on either side of the 
pressure minimum, the difference in pressure is twice the 
pressure at one of the coil locations. The pressure 
gradient is the pressure difference divided by the 
separation between the coils, and is given by the following 
formula: 

Vp = 2P /^x 
+ 

(4.6) 
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Gradient sensitivity is calculated directly by the 
■following -Formula, when the pressure gradient is sufficient 
to null the amplitude of the fundamental: 



M = 3.83rads/yp 
D 



(4.7) 



The units of gradient sensitivity are radians/ jj. Pa/cm. 

Gradient sensitivity can also be calculated indirectly 
by using single coil sensitivity values. Indirect gradient 
sensitivity can be expressed as follows: 



M 



ID 






(4.8) 



where W and W are the phase shifts in coil 1 and coil 2 
1 2 

respectively. Using the pressure value as expressed in 
equation (4.S) the phase shifts may be expressed as: 



<i = p * M 

1 + Cl 



(4.9) 



and 




P * M 

C2 



(4. 10) 
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where P and P are the pressure amplitudes at the sensor 
+ - 

coil locations and M and M are single coil sensitivities 

Cl C2 

obtained -from Table III. Inserting equations (4.9) and 

equation (4.10) into equation (4.8) and using the -Fact that 

P = -P when the sensor is at a pressure node, the indirect 
+ - 

gradient sensitivity can be expressed ast 

M = P CM + M 1/\7p 
ID ± Cl C2 

(4.11) 

Substituting equation (4.6) into equation (4.11) results in 
an indirect gradient sensitivity -Formula as -Follows: 

M = CM + M 3/C2/^x3 
ID Cl C2 

(4. 12) 

where is the separation between the two coils and this 

distance is assumed to be smaller than a wavelength. 

Table IV is a summary o-F direct and indirect gradient 
sensitivity measurements calculated -From dual coil and 
single coil data respectively. Measurement o-F direct 
gradient sensitivity was calculated at three di-F-Ferent 
calibrator resonances and di-F-Ferent depths associated with 
pressure minima (maximum pressure gradient). Indirect 
gradient sensitivity is also calculated at three di-F-Ferent 
calibrator resonances but uses the average sensitivity 
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obtained at a pressure maxima for the single coils divided 
by the desired distance between the coils (equation (4.12)). 

Error analysis for direct gradient sensitivity can be 
made by taking the square root of the sum at squares of the 
constituents of equation (4.6). This may be expressed as 
follows: ^ 

'dM\ 

(4.13) 



M 




where C is the experimental average speed of sound in the 
calibrator as shown in Table I. All other symbols are as 
explained in equation (4.3) and equation (4.5). 

Error analysis for indirect gradient sensitivity may be 
expressed as follows: 






M 



^ID 







where <M> is 



the 







(4. 14) 


for 


single 


coi 1 


the 


single 


coi 1 



sensitivities as presented in Table III< 



D. ANALYSIS 



A general analysis of single coil sensitivity compared 
to published data CRef.63 can be shown. From Table III 
typical sensitivities are on the order of approximately 
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-2 

10 iJLfad/ jJLPa per hydrophone for 10 meters of fiber. This 
results in a sensitivity of 10 flraLd/jJ. Pa/m. For an optical 
wavelength of 0.83 microns, the experimental value of 
A<PI4>Ap is approximately 130 x 10 pascal . For the 
same diameter fiber (0.4mm), the calculated value is greater 
than published values CRef. 51 by a factor of four to five. 
This variation could be explained in terms of the large 
variations in observed physical properties of Hytrel 
coatings resulting in a more sensitive fiber. 

As shown in Table III, the single coil sensitivities of 
coil 1 and coil 2 at the same depth appear to be closely 
matched within experimental error. The sensitivity of the 
individual coils are not matched over a depth range in the 
acoustic calibrator. The data indicates that the 
sensitivity of the single coil is a function of depth with 
possible near field effect occuring close to the source 
causing an anomaly. No reason could be found for this 
result but follow-on experiments will investigate this depth 
dependent sensitivity of single coils. 

Measured pressure gradient sensitivity and calculated 
pressure gradient sensitivity as shown in Table IV in the 
previous section are in agreement within experimental 
error. The measured pressure gradient sensitivity obtained 
from data at pressure minima is shown to be more constant 
with frequency and depth than the calculated pressure 
gradient sensitivity based on single coil measurements at 
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pressure maxima. Where the single coil data at Table III is 
a -function o-f depth the measured sensitivity o-F the dual 
coil setup appears to be relatively independent o-F depth 
with an average value o-F 0. 124+0. OIS ^rad/^ Pa/cm. 

Final analysis o-F the pressure gradient hydrophone was 
conducted by connecting the output o-F the photodiodes to the 
phase-lock loop locking circuit. The dual -Fiber coils with 
a LC-10 hydrophone in the center was lowered in the acoustic 
calibrator. The voltage from the locking circuit and the 
LC-10 hydrophone were recorded. Figures 4.6 and 4.7 are 
plots of the recorded data at 593 Hz and 752 Hz 
respectively. As shown in Figures 4.6 and 4.7 as LC-10 
hydrophone proceeds to a minimum gradient hydrophone voltage 
proceeds to a maximum and as LC-10 hydrophones proceeds to a 
maximum gradient hydrophone voltage proceed to a minimum. 
The data to produce Figures 4.6 and 4.7 is in good agreement 
with theory except for data recorded near the source or near 
the air — water boundary. 
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Figure 4.6 Demodulated Gradient Hydrophone 
Output at 593 Hz. 
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DISTANCE FROM TOP OF WATER COLUMN-CENTIMETERS 
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Figure 4.7 Demodulated Gradient Hydrophone 
Output at 752 Hz . 
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DISTANCE FROM TOP OF WATER COLUMN-CENTIMETERS 



V/. CONCLUSIONS AND RECOMMENDATIONS 



Since their invention in the late 1970's CRe-f. 3 and 41, 
many applications for fiber optic sensors have been 
suggested. In the area of low frequency (audio) sound 
transduction in the sea (hydrophones) , it was felt the most 
reasonable application of this class of sensor was as a 
dipole sensor. The potentially high sensitivity could be 
exploited since the output of a dipole sensor is the small 
difference of two relatively large signals. Secondly, an 
interf erometric (homodyne) detector is intrinsically 
differential since, in this case the interference is created 
by the diff eren ce of the optical path lengths in the two 
arms of the interf erometer . An additional benefit of this 
differential operation is the high degree of common mode 
noise rejection which maybe gained by having the "sensor'' 
and "reference" arm in the same environment until the very 
end of the sensor. 

The laboratory investigations outlined in this report 
have demonstrated the feasibility of this application. 
Because the fiber optic system is essentially a phase 
modulated rather than amplitude modulated system (equation 
2.22), it is difficult to make an absolute comparison 
between it and a piezoelectric pressure gradient sensor. 
Before the impedance converting, 20 dB gain, internal 
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amplifier, the USRD type G63 standard pressure gradient 

-9 

hydrophone CRef. 13D has a sensitivity of 7.5X10 V/^Pa/cm. 
The fiber optic system described here using piezoelectric 
phase shifter (Chapter III section A-4) and the feedback 
locking circuit (Chapter III section D-8) had an output of 
305+8 mV/radian (Chapter IV section A and Fig. 4.3). This 
combined with an average gradient sensitivity of 
0.124+0.015 ^ rad/// Pa/cm gives an overall sensitivity for 
the fiber system of 3.8+0.5X10 V/yL£ Pe/cm. By this method, 
the fiber system has a 14 dB advantage over the USRD type 
863. 

At this stage, further work in this area is justified. 
The directions that should be pursued are the ruggedizing of 
the sensor, and the development of a package which will 
incorporate the sensor and couplers in a single submersible 
unit, possibly to be followed by the integration of the 
laser and demodulators in the same unit. Alternative 
demodulation schemes which do not require a piezoelectric 
fiber stretcher to maintain "quadrature lock" should also be 
investigated. Finally, direct comparison with other 
directional dipole sensors such as those on DIFAR sonobuoys 
should be made under real sea conditions. 
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APPENDIX A 



DATA ACQUISITION PROGRAM 

This program was written to use a HP-85F computer to 
gather data more accurately and quickly than by manual 
means. This program was used after initial construction and 
prior to sensitivity experiments to obtain data -for 
computation of interf erometer characteristics (see Figure 
4.1). Chapter IV Section A gives a detailed explanation of 
how this program was utilized for computerized data 
acquisition. 
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10 OPTION BASE 1 

20 SHORT F(3), A(250) , B(250,3), 0(250,3), E(250) 

30 REM ** OBTAIN INITIAL VALUES FOR PARAMETERS ** 

40 DISP "ENTER FUND FREQ" 

50 INPUT F(l) 

60 DISP "ENTER SPECT ANAL START FREQ" 

70 INPUT FI 

80 DISP "ENTER START AMPLITUDE IN MVOLTS" 

90 INPUT A(l) 

100 DISP "ENTER END AMPLITUDE IN MVOLTS" 

110 INPUT VI 

120 DISP "ENTER INCREMENT VALUE IN MVOLTS" 

130 INPUT I 
140 F(2)»2*F(1) 

150 F(3)«3*F(1) 

160 REM ** PRESET AND INITIALIZE HP-3314 AND HP-3582 ** 

170 OUTPUT 707 ; "PR" 

180 OUTPUT 707 ; "FR" ,F < 1 ) , "HZAP, A ( 1 ) , "MVFUl " 

190 OUTPUT 711 ;"PRS" 

200 OUTPUT 711 j "AS8SL1MN1MD4SP10AD",F1,"SC1AV4NU4AM2" 

210 J=(V1-A(1) ) /I+l 

220 REM ** DATA ACQUISITION FOR SET DRIVE AMPLITUDE «« 

230 FOR R=1 TO J 
240 OUTPUT 711 } "RE" 

250 WAIT 11500 
260 ENTER 723 ; V 
270 E(R)+V 

280 REM ** STORE AMPLITUDE OF FUND AND 1ST TWO HARMONICS ** 

290 FOR N=1 TO 3 

300 M=(F(N)-F1)/1000*250 

310 OUTPUT 711 ;"MP",M 

320 OUTPUT 711 ; "LMK" 

330 ENTER 711 ; B (R,N) ,C (R,N) 

340 NEXT N 

350 D=R+1 

360 A(D)»A<R)-('I 

370 OUTPUT 707 ; " AP" , A (D) , "MV" 

380 NEXT R 

390 REM ** OUTPUT DATA TO PRINTER ** 

400 PRINTER IS 701,80 

410 PRINT "FUND FREQ (HZ) =",F(1) 

420 PRINT 'SOURCE VOLTS LCIO" 

430 PRINT "MILLIVOLTS (P-P) VOLTS(RMS) VFl (RMS) VF2(RMS) 
VF3(RMS) " 

440 FOR R=1 TO J 

450 IMAGE 50,12X,D.DDE,4X,D.DDE,4X,D.DDE,4X,D.DDE 
460 PRINT USING 400 ; A(R) ,E(R) ,B(R, 1) ,B(R,2) ,B(R,3) 
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470 NEXT R 

480 REM ** RESET INSTRUMENTATION TO PRESET CONDITION ♦♦ 

490 OUTPUT 707 ; "PR" 

500 OUTPUT 711 j"PRS" 

510 REM ** STORE DATA ON TAPE ** 

520 ASSIGN* 1 TO "HYDl" 

530 PRINT# 1,1 i AO 
540 PRINT# 1,2 } EO 
550 PRINT# 1,3 } B(,) 

560 ASSIGN# 1 TO * 

570 REM ** PLOT DATA ON COMPUTER PRINTER *♦ 

580 PRINTER IS 2 
590 FOR K=1 TO 3 
600 GCLEAR 

610 SCALE 0,10000,0,. 01 

620 XAXIS 0,1000 

630 YAXIS 0,.001 

640 PENUP 

650 FOR I»1 TO J 

660 PLOT A(I) ,B(I,K) 

670 NEXT I 
680 COPY 
690 NEXT K 
700 GCLEAR 

710 SCALE 0,10000,0,. 01 

720 XAXIS 0,1000 

730 YAXIS 0,.001 

740 FOR K=1 TO 3 

750 PENUP 

760 FOR 1=1 TO J 

770 PLOT Ad) ,B(I,K) 

780 NEXT I 
790 NEXT K 
800 COPY 
810 DISP "END" 

820 END 
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APPENDIX B 



SAMPLE RAW DATA 



PZT #1 

Laser Diode 35 mi 11 lamps 
Current 



Fundamental 

■frequency (Hz) 372 



Source 


Fundamental 


1st Harmonic 


2nd Harmon ii 


Voltage 


Voltage 


Voltage 


Voltage 


(volts) 


(millivolts) 


(millivolts) 


(mi 1 1 i volts 


(rms) 


(rms) 


(rms) 


(rms) 


0.035 


0.783 


0.030 


0.0019 


0.053 


1.180 


0.066 


0.0039 


0.071 


1.570 


0. 105 


0. 0048 


0.071 


1.980 


0. 182 


0.0106 


0.089 


2.320 


0.239 


0.0164 


0. 106 


2.700 


0.349 


0.0241 


0. 124 


3.050 


0.456 


0.0347 


0. 142 


3.430 


0.544 


0.0473 


0. 159 


3.800 


0.642 


0.0647 


0. 177 


4.070 


0.812 


0.0830 


0. 195 


4.380 


0.916 


0.1150 


0.230 


4.820 


1.130 


0. 1460 


0.247 


5.010 


1.270 


0. 1720 


0.247 


5.340 


1.510 


0.2060 


0.265 


5.590 


1.610 


0.2540 


0.283 


5.790 


1.750 


0.2880 


0.300 


6.110 


2.010 


0.3280 


0.318 


6.310 


2. 150 


0.4070 


0.336 


6.460 


2.370 


0.4520 


0.353 


6.850 


2.440 


0.5200 


0. 373 


6.870 


2.460 


0.5770 


0.390 


7. 130 


2.960 


0.6490 


0.409 


7.240 


2.970 


0.7160 


0.426 


7.330 


3.320 


0.8160 


0.444 


7.370 


3.420 


0.8880 


0.461 


7.550 


3.630 


0.9650 


0.478 


7.700 


3.830 


1 . 0500 


0.497 


7.750 


3.990 


1.1300 
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0 . S 13 


7.650 


4.070 


1 . 2600 


0.532 


7.750 


4.320 


1.3100 


0.549 


7.890 


4.250 


1 . 4400 


0.567 


7.460 


3.990 


1 . 5000 


0.584 


7.730 


4.690 


1 . 6700 


0.603 


7.770 


4.360 


1 . 7400 


0.620 


7.650 


5.030 


1.8100 


0.638 


7.450 


4.720 


1 . 9400 


0.655 


7.460 


5. 190 


2.0200 


0.672 


7.350 


4.910 


2.2200 


0.691 


7.280 


5.080 


2.2900 


0.707 


7.070 


5.450 


2. 4000 


0.726 


6.930 


5.460 


2.5400 


0.743 


6.850 


5.650 


2.6700 


0.761 


6.590 


5.730 


2.8100 


0.778 


6.500 


5.820 


2.9500 


0.797 


6.200 


5.970 


3.0400 


0.814 


6.050 


6.050 


3. 1800 


0.832 


5.760 


6 . 100 


3.2900 


0.849 


5.550 


6.320 


3.3600 


0.868 


5.290 


6.230 


3.4900 


0.885 


5.070 


6.370 


3.6100 


0.902 


4.770 


6.500 


3.6900 


0.920 


4.490 


6.480 


3.8500 


0.937 


4.200 


6.540 


4.0000 


0.956 


3.870 


6.540 


4.0600 


0.973 


3.590 


6.480 


4. 1900 


0.991 


3.270 


6.480 


4.2200 


1.010 


3.000 


6.560 


4.4500 


1.030 


2.710 


6.500 


4. 5600 


1.040 


2.360 


5.910 


4.6200 


1.060 


2.040 


6.210 


4.7300 


1.080 


1.750 


6 . 170 


4.7200 


1.100 


1.480 


6 . 120 


4.8400 


1.110 


1.170 


6.070 


4.8200 


1.130 


0.858 


5.550 


5.0100 


1.150 


0.617 


5.630 


5. 1200 


1.170 


0.346 


5.460 


5. 0600 


1.190 


0.339 


5.320 


5. 1600 


1.200 


0.611 


5 . 150 


5.2500 


1.220 


0.841 


4.840 


5.2900 


1.240 


1.080 


4.710 


5.2800 


1.260 


1.390 


4.570 


5.2900 


1 . 270 


1.610 


4.280 


5.2800 


1.290 


1.920 


4.050 


5.2900 


1.310 


2 . 170 


3.630 


5.3600 


1.330 


2.370 


3.560 


5.4000 


1.340 


2.600 


3.370 


5.3400 


1.360 


2.830 


3.070 


5.2500 


1.380 


2.980 


2.840 


5.2500 


1.400 


3.260 


2.660 


5. 1900 



100 



1.420 


3.340 


2.410 


5. 1400 


1.430 


3.540 


2.210 
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